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HYCOM  2.2  (I) 


First  public  release  of  HYCOM  2.2 
o  Scheduled  for  December,  2004 
Maintain  all  features  of  HYCOM  2.1 

o  Orthogonal  curvilinear  grids 

o  Can  emulate  Z  or  Sigma  or  Sigma-Z  models 

o  Explicit  support  for  1-D  and  2-D  domains 

o  KPP  or  Kraus-Turner  or 
Mellor-Yamada  2.5  or  Price-Weller-Pinkel 

o  Rivers  as  bogused  surface  precipitation 

o  Multiple  tracers 

o  Off-line  one-way  nesting 

o  Scalability  via  OpenMP  or  MPI  or  both 

*  Bit-for-bit  multi-cpu  reproducibility 

New  diagnostics  within  HYCOM 

o  Time-averaged  fields  (in  archive  files) 
o  Drifters 


HYCOM  2.2  (II) 

Alternative  scalar  advection  techniques 

o  Donor  Cell,  FCT  (2nd  and  4th  order),  MPDATA 

Vertical  coordinate  changes 

o  Vertical  remapping  uses  PLM  for  fixed 
coordinate  layers 

o  Thin  deep  iso-pycnal  layers 

o  Spatially  varying  iso-pycnal  layer  target  densities 

o  Stability  from  locally  referenced  potential  density 

Atmospheric  forcing  changes 

o  Option  to  input  ustar  fields 
o  Option  to  relax  to  observed  SST  fields 
o  Improved  COARE  3.0  bulk  exchange  coefficients 
o  Black-body  correction  to  longwave  flux 

Mixed  layer  changes 

o  GISS  mixed  layer  model 
o  KPP  bottom  boundary  layer 
o  KPP  tuning 

o  Latitudinally  dependent  background  diffusion 


HYCOM  2.2  (III) 

Improved  support  for  rivers 

o  Still  bogused  surface  precipitation 
o  Better  control  of  low  salinity  profiles 
o  Option  for  mass  (vs  salinity)  flux 

Nesting  no  longer  requires  co-located  grids 

o  General  archive  to  archive  horizontal 
interpolation 

Hybrid  to  fixed  vertical  grid  remapper 

o  Allows  fixed-coordinate  nests  inside  hybrid 
coordinate  outer  domains 

*  HYCOM  to  (fixed-grid)  HYCOM 

*  HYCOM  to  NCOM 

Diagnostic  fields  to  netCDF  and  other  file  formats 

o  All  x-y  “hycomproc”  fields 

*  Layer  space 

*  Velocity  interpolated  to  the  p-grid 

o  All  3-D  archive  fields  interpolated  to  z-space 

*  On  p-grid,  or 

*  Sampled  along  arbitrary  tracks 
o  Forcing  input  fields 


HYCOM  CURVI-LINEAR  GRIDS  and  NetCDF 

Most  basin-scale  cases  use  a  Mercator  grid 

o  l-D  latitude  and  longitude  axes 
o  Handled  well  by  many  netCDF  packages 

Global  HYCOM’s  Arctic  patch  grid  is  curvi-linear 

HYCOM  netCDF  use  the  CF-1 .0  conventions,  which 
support  curvi-linear  grids 

o  If  latitude  and  longitude  are  2-D  grids 

*  1-D  axes  are  array  indexes 

*  Longitude  and  latitude  arrays  are  also  in  the 
file  and  identified  as  alternative  coordinates 

Most  netCDF  packages  are  not  CF-1 .0  aware 

o  Can  plot  in  “logical”  (array)  space 

o  Interpolate  to  a  1-D  latitude  and  longitude 
grid  off-line 

*  General  archive  to  archive  horizontal 
interpolation 

Archive  to  archive  remapper  can  also  be  used  for 
standard  (non-native)  grids 


o  Mersea  grid  is  uniform  1/8° 


GoM  NESTED  TEST  DOMAIN 
Same  resolution  nesting  unexpectedly  useful 

o  No  need  to  rerun  large  domain 
o  Change  atmospheric  forcing  (e.g.  use  MM5) 
o  Change  vertical  structure 
o  Tracer  studies  (e.g.  add  biology) 

1/12°:  Gulf  of  Mexico  inside  Atlantic 

o  Change  from  20m  to  5m  coastline 
o  Run  for  Aug  1999  to  equilibrate 
o  Run  Sep-Nov  as  standard  test  case 

Used  to  test  advection  schemes 


ATLANTIC  vs  GOM  NEST  (MPDATA) 


MPDATA  VS  LEAPFROG-FCT  (SSH) 


MPDATA  VS  LEAPFROG-FCT  (SSS) 


layer=01  salinity  Dec  01, 1999  00Z  [03.1  H] 


98W  96W  94W  92W  90W  88W  86W  84W  82W  80W  78W 


CANDIDATE  FEATURES  FOR  HYCOM  2.3 
Stable-code  vs  new  features 

o  Released  code-base  has  to  be  tested  and  stable 
o  New  features  can  be  a  significant  improvement 
o  Will  add  interim  releases  to  web  page 
*  Features  may  be  removed  in  next  released  code 

Fully  region-independent 

o  Compile  once,  run  on  any  region  and  any 
number  of  processors 

o  Needed  for  full  ESMF  compliance 
Improve  split-explicit  time  scheme 
Tidal  forcing 
Diurnal  heat  flux  cycle 
Equation  of  state  that  is  quadratic  in  salinity 
Even  better  support  for  rivers 
Wind  drag  coefficient  based  on  model  SST 
Initial  support  for  ESMF 


HYCOM  AND  ESMF 

Earth  System  Modeling  Framework 
http://www.esmf.ucar.edu/ 

o  Superstructure  couples  components 

*  Air/Ocean/Ice/Land 

*  Asynchronous  I/O  component 

•  Not  yet  available  via  ESMF 

o  Infrastructure  provides  data  structures 
and  utilities  for  building  scalable  models 

Add  a  superstructure  “cap”  to  HYCOM 

o  Simplifies  coupled  systems 

*  HYCOM  coupled  to  LANL  CICE  sea-ice 

*  Convert  atmospheric  field  processing  and  the 
energy-loan  ice  model  into  ESMF  components 

o  Use  ESMF  for  I/O 

This  initial  ESMF  support  will  probably  be  optional 

ESMF  may  be  required  to  run  HYCOM  at  some  point 

o  Harder  to  get  started  with  HYCOM 
o  Will  provide  many  new  capabilities 


HYCOM  AND  HOME 

Hybrid  Ocean  Modeling  Environment  (HOME) 

o  Not  one  model,  but  an  environment 

o  Unify  existing  isopycnal/hybrid  ocean 
models  into  a  single  code  base 

o  Still  an  unfunded  proposal 

There  will  be  a  migration  path  from  HYCOM  to  HOME 

o  Re-implement  HYCOM  in  HOME 
o  HYCOM  with  ESMF  will  simplify  the  migration 

HOME  “best  practices”  studies  may  find  better 
alternatives  to  HYCOM  algorithms 

o  Exact  mass  conservation 
o  Better  free  surface  formulation 
o  Improved  time  stepping 

Some  of  these  may  be  back  ported  to  HYCOM 

At  some  point  “HYCOM  in  HOME”  will  become 
the  only  supported  HYCOM 

o  Might  be  very  different  to  HYCOM  2.X 
o  Might  not  even  be  called  HYCOM 


